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Abstract. - Most real systems are growing. In order to model the evolution of real systems, 
many growing network models have been proposed to reproduce some specific topology prop- 
erties. As the structure strongly influences the network function, designing the function- aimed 
growing strategy is also a significant task with many potential applications. In this letter, we 
focus on synchronization in the growing networks. In order to enhance the synchronizability 
during the network evolution, we propose the Spectral-Based Growing (SBG) strategy. Based 
on the linear stability analysis of synchronization, we show that our growing mechanism yields 
better synchronizability than the existing topology-aimed growing strategies in both artificial and 
real-world networks. We also observe that there is an optimal degree of new added nodes, which 
means adding nodes with neither too large nor too low degree could enhance the synchronizability. 
Furthermore, some topology measurements are considered in the resultant networks. The results 
show that the degree, node betweenness centrality from SBG strategy are more homogenous than 
those from other growing strategies. Our work highlights the importance of the function- aimed 
growth of the networks and deepens our understanding of it. 



Introduction. — Networks, despite their simplicity, 
represent the interaction structure among components in 
a wide range of real complex systems, from social rela- 
tionships among individuals, to interactions of proteins in 
biological systems, and even the interdependence of func- 
tion calls in large software projects. The network con- 
cept has been developed as an important tool for analyz- 
ing the relationship of structure and function for many 
complex systems in the last decades [IH5]. An interesting 
phenomenon observed in complex networks is synchroniza- 
tion [I]. As a kind of collective behavior, synchronization 
is often encountered in living systems, such as circadian 
rhythm, phase locking respiration with mechanical ven- 
tilator, phase locking of chicken embryo heart cells with 
external stimuli [5]. The early works on synchronization 
were concerned with only a small number of coupled os- 
cillators. However, as the synchronization in many real- 
world systems is based on a large number of dynamical 
units interacting with a complex coupling structure, the 
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research of synchronization on complex networks has at- 
tracted much attention recent years. 

Previous studies have given us an overview over the syn- 
chronizability of some well-known network models. In gen- 
eral, random networks have better synchronizability than 
regular networks, and Strogatz- Watts networks have bet- 
ter synchronizability than scale- free networks [7r 9^ . Actu- 
ally, there are many factors that affect the network syn- 
chronizability. For example, average shortest distance is 
an important factor. Besides, the heterogeneity of the net- 
work is also one of the most influential factors determining 
the synchronizability [TU]. Generally speaking, the less 
heterogeneous the network is, the better its synchroniz- 
ability will be. With these understandings, many meth- 
ods have been proposed to enhance synchronizability in 
complex networks. There are several main directions in- 
cluding designing strategies for coupling strength [TTHTB] . 
modifying network structures |17fl20j . flipping the direc- 
tionality pTH25] and so on. Each group of methods has 
its corresponding application field since the properties of 
real systems may form different operation constraints. 
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Most of the former works on synchronization focus on 
given size networks. However, most real- world networks 
are growing, which is evidenced by numerous systems in- 
cluding man-made and natural ones. Most of the previ- 
ous growing mechanisms aim at reproducing some spe- 
cific topology properties. For instance, in the pioneering 
work of Barabasi- Albert (BA) network [3], the preferen- 
tial attachment growing mechanism reproduces the power- 
law degree distribution observed in many realistic complex 
systems. However, the growing mechanisms for network 
functions also need to be explored. For example, in some 
technological networks such as electric power grids, the 
neural networks and even social networks, the introduc- 
tion of new nodes are supposed to enhance the network 
functions [23J[2S]. In these cases, the topology properties 
of the network should be by-products during the improve- 
ment of the functions [5^ . Therefore, it is interesting to 
study the strategy to enhance functions in the growing 
networks. 

In this letter, employing synchronization as the typical 
function, we propose the Spectral-Based Growing (SBG) 
strategy. We compare the new growing mechanisms with 
some existing topology-aimed growing mechanisms includ- 
ing the Preferential Attachment (PA), Reversed Preferen- 
tial Attachment (RPA), Random Attachment (RA) mech- 
anisms. We find that the SBG strategy yields a better syn- 
chronizability in both artificial and real-world networks. 
Moreover, detailed study on some topology measurements 
of the resultant networks shows that the topology mea- 
surement such as the degree and node betweenness cen- 
trality in SBG networks are more homogeneous than those 
in the networks from other growing mechanisms. Besides, 
the SBG networks enjoy smaller clustering coefficient and 
average shortest path length. Interestingly, we observe 
that the synchronizability enhancement is strongly related 
to the degree of new added nodes. The synchronizability 
would be weaken if the degree of the new nodes are too 
large or too low. Actually, there is an optimal degree 
for each given system and the optimal degree is approx- 
imately proportional to the average degree of the origi- 
nal networks. Finally, in order to overcome the drawback 
of trapping in local optimality of SBG method when a 
large number of nodes are added, we tested the simulated- 
annealing-based SBG (SSBG) strategy. With the simu- 
lated annealing process, the SSBG strategy can jump out 
of the local optimality and further improve the synchro- 
nizability when a large number of nodes are added to the 
networks. Our work deepens our understanding of the 
function-aimed growth of networks and may have wide 
potential applications. 

The rest of the paper is organized as follows. In Sec. 
II, a new growing mechanism called the SBG strategy is 
proposed for enhancing the synchronizability of growing 
networks. In Sec. HI, we will detailedly investigate the 
topology of the SBG networks and the optimal degree for 
new nodes. The application to real-world networks is in 
Sec. IV. In Sec. V, the SBG strategy with simulated an- 



nealing process is also tested. Finally, concluding remarks 
are presented in section VI. 

Spectral-Based Growing (SBG) Strategy. 

Let's consider a system with N identical oscillators sym- 
metrically coupled through a network. The equations of 
motion for the oscillator state vector x at each node i are 

N 

±1 = F{xi) ~a^L.ijH{xj), (1) 

where F{x) determines the uncoupled oscillator dynam- 
ics of each node, H{x) specifies the coupling of the vec- 
tor fields, a is the coupling strength, and L is the graph 
Laplacian. Specifically, L is defined as Lij = —1 if an edge 
exists between node i and j, and La = ki, where ki is the 
degree of node i. Solutions of this equation are defined 
to be synchronized if Xi = Xj for all nodes i and j in the 
network. From ref. [9j, for many oscillatory systems the 
master equation is negative only in a single interval [ai, 02] 
determined by F and H. This implies that the network is 
synchronizable only when the eigenratio Xn/^2 < dijotx. 
One can conclude that the synchronizability is related to 
the eigenratio r ~ A7V/A2, where A^r and A2 are the largest 
and second-smallest eigenvalues of the graph Laplacian, 
respectively. The smaller the eigenratio r of a network 
is, the better its synchronizability will be. Our aim is 
to enhance the synchronizability when the network grows, 
which is equivalent to minimizing the eigenratio r. 

For a symmetric matrix L, the eigenvalues of this ma- 
trix and the corresponding eigenvectors are denoted as A 
and V respectively. When an edge is added to or removed 
from the network, we denote AL, AA and At; as the cor- 
responding perturbation of the Laplacian matrix L, its 
eigenvalues and its eigenvectors. We have the following 
equations: 

(L + AL)(ii+ At;) = (A + AA)(i;-f At;). (2) 

Left multiplying (2) by iF and neglecting second order 



terms v^ AL5v and v^ A\5v, we obtain AA, 



^ AL^ 



If 



an edge between node m and n is added to the network, 
^Lmm = AL„„ = 1 and ALm„ = AL„„i = -1. Then 



Ai increases by {vi,r. 



if the edge is added to the 



network, where Vi is the eigenvector with unit norm cor- 
responding to \i . This gives the first-order approximation 
of the increase in A^ . 

When an edge is going to be removed from the net- 

\ aT"^^^ is satisfied, we choose 
A2— AA2.6 ' 



work, if ^"'^^"■° < ^ 



A2-AA2 

edge a instead of edge b since the eigenratio after re- 
moving edge a is smaller. Neglecting the second-order 
terms, we obtain that an existing edge a which minimizes 
AArAA2,a — A2AAAr_a should be removed and a nonexist- 
ing edge a which maximizes AArAA2,a — A2AAAr,a should 
be added. Accordingly, we propose the Spectral-Based 
Growing Strategy (SBG) as follows: 

1. A new node i with k links is added, and these k links 
are connected randomly to the original network. 
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2. For each existing edge connecting node i and node 
j of the network, the quantity C — Ajv(w2,j — "^2,^)^ ^ 
^2{vN.i — VN.jY is calculated. The edge with minimum C 
is cut. 

3. For each pair of unconnected nodes i and I of the 
network, the quantity i? = \N(v2,i—V2,iY—^2(vN,i—VN,iY 
is calculated. Then an edge between a pair of nodes with 
maximum E is created. 

4. If the eigenratio Xn/M of the new network becomes 
smaller after steps 2 and 3, the rewiring is accepted, oth- 
erwise the process is rejected. 

5. Once step 4 is rejected, all the links for this new node 
i are fixed. 

We compare the SBG strategy with some topology- 
aimed growing strategies including the preferential attach- 
ment (PA), reversed preferential attachment (RPA) and 
Random attachment (RA) strategies. In PA, the proba- 
bility of a node i in the network connected by the new 
node is proportional to ki, where ki is the degree of node 
i. This growing mechanism yields the power-law degree 
distribution observed in many realistic systems. In RPA, 
the probability of a node i in the network connected by 
the new node is proportional to k^ . Consequently, the 
node degree will be homogenous. In RA, the new nodes 
randomly connect to existing nodes in the original net- 
work. 

In the following simulation, we mainly consider four dif- 
ferent kinds of network models: Barabasi- Albert (BA) net- 
work with fc = 10 [3 ; Watts-Strogatz (WS) small-world 
network with p = OM and A: = 20 [T; Erdos-Renyi (ER) 
random graph with A: = 15 [27;; and Girvan-Newman(GN) 
benchmark with A;,;„ — 12 and k = 16 [28| . The initial 
number of nodes in BA, SW and ER modeled networks is 
300, and the size of GN benchmark is 128. When a new 
node is introduced, the degree of it is equal to the average 
degree of the original network. Fig. [l]shows the eigenratio 
r as a function of the number of nodes added. Obviously, 
the eigenratio in SBG strategy decreases fastest in all these 
four network models. The eigenratio in RPA strategy per- 
forms second best because it yields a homogeneous degree 
distribution which is favorable for synchronization. As 
a widely used growing strategy, the PA strategy cannot 
lead to a good synchronizability since heterogenous de- 
gree distribution is formed by the PA strategy. Our result 
indicates that the strategy aiming at reproducing topol- 
ogy properties may not guarantee a high-performance of 
some typical function. 

Topology of the SBG networks and the optimal 
degree for new nodes. — In this section, we will de- 
tailedly discuss topology properties of the SBG networks 
and the effect of new nodes' degree on the resultant syn- 
chronizability. 

First, we investigate the topology properties of the net- 
works after the growth including degree, clustering coeffi- 
cient, average shortest path length, node betweenness cen- 
trality. Starting from the initial networks with 300 nodes. 
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Fig. 1: (color on-line) Four growing mechanisms in different 
network models including (a) ER random graph with k = 15, 
(b) GN benchmark with kin ~ 12, (c) WS small-world network 
with p = 0.01 and k = 20, (d)BA network with k = 10. The 
initial network size in these four models is 300. The results are 
obtained by averaging over 20 independent realizations. 



20 new nodes are added and the degree of them is equal 
to the average degree of the original network. We com- 
pare SBG networks with the initial, PA, RPA, and RA 
networks and the results are shown in table [TJ From the 
results, we observe that the standard deviation of node 
degYee{VARk), the average clustering coefficient(< C >), 
the average shortest path length(< d >), the standard 
deviation of node betweenness centrality(V'Ai?h) and the 
maximum node betweenness centrality(Mj4X5) from the 
SBG strategy are generally smaller than those from the 
topology-aimed growing strategies. It has been pointed 
out that networks with optimal synchronizability should 
be a class of entangled networks with homogenous topol- 
ogy properties, such as degree, node betweenness central- 
ity [17]. From the results in table [l] we can see that the 
SBG networks are closest to the optimal state among all 
these growing strategies. 

Second, we find that the synchronizability enhancement 
is strongly related to the degree of new added nodes. If 
the degree of new added node is too small or too large, 
the synchronizability cannot be enhanced. Consider a ER 
random graph with the size of the network iV = 100 and 
average degree k = 5, the eigenratio of it is r = 17.96. If 
a new node is introduced with node degree A: = 1 , the ex- 
haustive search method finds that the smallest eigenratio 
is r = 18.01 which is even larger than that in the original 
network. If a new node is added with degree k = 100, the 
eigenratio is r = 55.365 which is also larger than that in 
the original network. That is to say, if the degree of new 
added nodes is too small (e.g. k — 1) or too large (e.g. 
k = N), the synchronizability would be weaken. There- 
fore, there should be an optimal degree for the new added 
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Table 1: The topology properties in the resultant networks from different growing mechanisms. 



Growing 


VARk 


<C> 


< d > 


MAXb 


VARb 


Original 


Mechanism 












Network 


SBG 


3.161 


0.051 


2.401 


1112.7 


184.04 




RPA 


3.568 


0.055 


2.403 


1266.3 


207.23 




PA 


3.910 


0.055 


2.402 


1315.6 


222.13 


ER 


RA 


3.736 


0.055 


2.403 


1301.7 


210.45 




SBG 


7.710 


0.076 


2.611 


5668 


887.17 




RPA 


7.600 


0.075 


2.621 


5912 


882.30 




PA 


7.950 


0.072 


2.601 


6434 


944.26 


BA 


RA 


7.870 


0.078 


2.612 


6579 


921.53 




SBG 


2.441 


0.180 


2.010 


274.8 


42.51 




RPA 


3.060 


0.190 


2.021 


321.9 


61.48 




PA 


3.860 


0.200 


2.033 


403.7 


75.55 


GN 


RA 


3.310 


0.182 


2.014 


405.8 


67.10 




SBG 


1.550 


0.350 


2.680 


2091.1 


273.3 




RPA 


1.780 


0.358 


2.690 


2000.0 


298.1 




PA 


2.017 


0.364 


2.700 


2392.8 


322.1 


SW 


RA 


1.847 


0.361 


2.700 


2338.5 


310.4 





nodes. We consider two situations in which one node and 
ten nodes are added to ER random networks, respectively. 
Fig. [JJa) and (b) shows the eigenratio as the function of 
the node degree of new added nodes. The results indicate 
that there is an optimal degree which minimizes the syn- 
chronizability of the resultant networks. In Fig. [2jc) and 
(d), we can see that the optimal degree of new added nodes 
is approximately proportional to the average degree of the 
original ER networks when using SBG method. Moreover, 
we remark that the similar relation is observed from the 
simulations on SW and BA models. 

The application to real- world networks. In this 
section, we will apply the SBG strategy to some real-world 
networks. We mainly consider four different kinds of net- 
works including neural, metabolic, power grid and social 
networks. Actually, synchronization has been shown to be 
of special relevance in these real systems [6] . Specifically, 
large-scale synchronization of oscillatory neural activity is 
believed to play a crucial role in the information and cog- 
nitive processing. In the metabolic networks, synchroniza- 
tion is important for chemical reactions of metabolism as 
well as the regulatory interactions that guide these reac- 
tions. Synchronization of the power grids is understood as 
every station and piece of equipment running on the same 
clock, which is crucial for its proper operation. Cascading 
failures related to de-synchronization can lead to massive 
power blackouts. In social systems, synchronization is of 
widely interested and related to many collective processes, 
such as opinion formation. 

Here, we employed C.elegans neural and metabolic net- 
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Fig. 2: (color on-line) The eigenratio as a function of the degree 
of new added nodes, (a) One node is added to the network, 
(b) Ten nodes are added to the network, (c) and (d) show the 
optimal degree in the SBG method as a function of the average 
degree of the original networks. The initial networks are ER 
models with A^ = 100. In (a) and (b), the degree of initial 
networks is fc = 10. The results are obtained by averaging over 
20 independent realizations. 



works, power grid of western US and the collaboration 
network among scientists working at the Santa Fe Insti- 
tute (SFI network) ^?|. We adopted the largest connected 
component of these networks and consider all of them as 
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Fig. 3: (color on-line) The eigenratio as a function of the de- 
gree of new added nodes when SBG strategy is applied to 
real-world networks. The real-world networks includes (a) 
C.elegans metabolic network; (b) power grid; (c) SFI network; 
(d) C.elegans neural network. New nodes are introduced with 
degree equal to the average degree of the original networks. 
The results are obtained by averaging over 20 independent re- 
alizations. 



undirected ones. The neural network of C. elegans con- 
tains 302 neurons and 2, 359 links with nodes represent- 
ing neurons and links representing synaptic connections. 
In the metabolic network, nodes represent the chemicals 
produced and consumed by the reactions. If a pair of 
chemicals are in the same reaction, then there is a link 
between them. C.elegans metabolic network contains 453 
vertices and 4, 596 edges. The power grid is the network of 
high- voltage transmission lines that provide long-distance 
transport of electric power in western US. The power grid 
network here is with 4, 994 nodes and 6, 594 links. The 
scientist collaboration network represents the collabora- 
tion relations among scientists working at the Santa Fe 
Institute. Edges are placed between scientists who have 
published at least one paper together. In SFI network, 
there are 118 nodes and 200 links. 

When new chemicals in metabolic networks, new neu- 
rons in neuron networks, new power plants in power grid 
or new comers in social systems are introduced, how to en- 
hance the network synchronizability is an important prob- 
lem which is meaningful in both theoretical and practical 
senses. Based on these real-world networks, twenty new 
nodes are introduced with degree equal to the average de- 
gree of the original networks. Fig. [3] shows the eigenra- 
tio r as a function of the number of new nodes. As ex- 
pected, we find that the SBG method has advantage in 
enhancing synchronizability compared to these topology- 
aimed methods in all the real- world networks. Moreover, 
the errorbar of the synchronizability from the SBG strat- 
egy is smaller than that from other methods. That is to 



say, our method is more effective and reliable. Just like 
our simulation on artificial networks, we also observe that 
RPA strategy leads to a good synchronizability among all 
the topology-aimed strategy. This result proves again the 
importance of homogeneous topology for synchronization 
when the network is growing. 

The SBG strategy with simulated annealing pro- 
cess. — When just a few nodes are added to a given 
network, SBG is an effective growing mechanism to en- 
hance the network synchronizability. However, the SBG 
strategy can be sometimes trapped in the local maxi- 
mum, especially when many nodes are added. For ex- 
ample, when the initial networks is with only 20 nodes, 
after adding 500 nodes the synchronizability from SBG 
strategy is similar to that from RPA strategy. In this sec- 
tion, we tested an extended SBG strategy which is called 
simulated-annealing-based SBG (SSBG) strategy. With 
the simulated annealing process, the SSBG strategy can 
jump out of the local maximum and further improve the 
SBG method when a large number of nodes are added 
to the networks. The algorithm consists of the following 
steps: 

1. A new node i with k links is added, and these k links 
are connected randomly to the original network. 

2. For each existing edge connecting node i and node 
j of the network, the quantity C — XN{v2,i — V2.jY — 
X2{vN,i — VN,if' is calculated. The probability for remov- 
ing an edge is proportional to ec. 

3. For each pair of unconnected nodes i and I of the 
network, the quantity -E = \N{v2,i—V2.iY ~X2{vn ,i—VN ,iY 
is calculated. The probability for creating an edge between 
the unconnected nodes is proportional to e^. 

4. If the eigenratio \n 1^2 of the new network be- 
comes smaller after steps 2 and 3, the rewiring is ac- 
cepted, otherwise the rewiring is accepted with the prob- 
ability of e("((^"/^2)N..„-(Ajv/A2)oM)/Tffi?0. Here, we fixed 
THR = 3. 

5. The process is stopped after 150 times of iterations. 
In our simulation, we start from a small ER random 

graph with N = 20, k = 10. Five hundred new nodes are 
introduced one by one. The degree of new added nodes 
is equal to the average degree of the original networks. 
Consequently, the resulting network is independent of the 
initial network since the size of the initial network is suf- 
ficiently small. Fig. |4] shows eigenratio r as a function of 
the number of new nodes added. The results indicates 
that SSBG method indeed yields a network with better 
synchronizability than the topology-aimed methods. 

Finally, from the computational point of view, the 
SSBG method is with higher computational complexity 
than the original SBG method. However, compared with 
the ordinary simulated-annealing method which directly 
enhances synchronizability by randomly moving one link 
in each iteration, SSBG method can significantly reduce 
the computation cost since steps 2 and 3 can speed up the 
convergence of the whole process. 
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Fig. 4: (color on-line) The eigenratio Xn/^2 as a function of 
the number of new nodes added. The initial network is an ER 
random graph with Ai" = 20, fc = 10. New nodes are intro- 
duced with degree equal to the average degree of the original 
networks. 



Conclusion. — In this letter, aiming at improving the 
synchronizability in the growing networks, we propose the 
Spectral-Based Growing (SBG) strategy. We compare the 
SBG strategy with some existing topology-aimed grow- 
ing ones including the Preferential Attachment (PA), Re- 
versed Preferential Attachment (RPA), Random Attach- 
ment (RA) mechanisms. The results show that the SBG 
strategy yields a better synchronizability in both artificial 
and real- world networks. 

Moreover, some topology measurements are considered 
in the resultant networks. It shows that the topology mea- 
surement such as degree and node betweenness centrality 
in SBG networks are more homogeneous than those in the 
networks from the other growing mechanisms. Interest- 
ingly, we find that the synchronizability enhancement is 
strongly related to the degree of new added nodes. The 
synchronizability would be weaken if the degree of the new 
nodes are too large or too low. Actually, there is an opti- 
mal degree for each given system and the optimal degree 
is approximately proportional to the average degree of the 
original networks. 

Finally, in order to overcome the drawback of trapping 
in local maximum of SBG method when a large number of 
nodes are added, we tested the simulated-annealing-based 
SBG (SSBG) strategy. With the simulated annealing pro- 
cess, the SSBG strategy can jump out of the local maxi- 
mum and further improve the SBG method when a large 
number of nodes are added to the networks. All these 
findings highlight the importance of the functional growth 
of networks and deepen our understanding of it. 
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